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Single photons are generated from an ensemble of cold Cs atoms via the protocol of Duan et al. 
[Nature 414, 413 (2001)]. Conditioned upon an initial detection from field 1 at 852 nm, a photon 
in field 2 at 894 nm is produced in a controlled fashion from excitation stored within the atomic 
ensemble. The single-quantum character of the field 2 is demonstrated by the violation of a Cauchy- 
Schwarz inequality, namely m(12, I2II1) = 0.24±0.05 ^ 1, where w(l2, 12II1) describes detection of 
two events (I2, 12) conditioned upon an initial detection li, with w ^ for single photons. 



A critical capability for quantum computation and 
communication is the controlled generation of single- 
photon pulses into well-defined spatial and temporal 
modes of the electromagnetic field. Indeed, early work on 
the realization of quantum computation utilized single- 
photon pulses as quantum bits {fiying quhits), with non- 
linear interactions mediated by an appropriate atomic 
medium 0, More recently, a scheme for quantum 
computation by way of linear optics and photoelectric de- 
tection has been developed that again relies upon single- 
photon pulses as qubits j^]- Protocols for the implemen- 
tation of quantum cryptography Q and of distributed 
quantum networks also rely on this capability j5|, ||6j|, as 
do some models for scalable quantum computation Q . 

Efforts to generate single-photon wavepackets can be 
broadly divided into techniques that provide photons "on 
demand" (e.g., quantum dots coupled to microcavities 

0. 0.0) and those that produce photons as a result of 
conditional measurement on a correlated quantum sys- 
tem. For conditional generation, the detection of one 
photon from a correlated pair results in a one-photon 
state for the second photon, as was first achieved us- 
ing "twin" photons from atomic cascades [ill fl^ and 
parametric down conversion 13] , with many modern ex- 
tensions [11 El El Eli. Within the context of the col- 
lective enhancement of atom-photon interactions in op- 
tically thick atomic samples |18L Il9| , a remarkable pro- 
tocol for scalable quantum networks suggests a new 
avenue for producing single photons via conditional mea- 
surement. 

Inspired by the protocol of Ref . [^ , in this Letter we 
report a significant advance in the creation of single pho- 
tons for diverse applications in quantum information sci- 
ence, namely the generation and storage of single quanta 
from an atomic ensemble. As illustrated in Figure 1, an 
initial write pulse of (classical) light creates a state of col- 
lective excitation in an ensemble of cold atoms as deter- 
mined by photoelectric detection for the generated field 

1. Although this first step is probabilistic, its success 
heralds the preparation of one excitation stored within 
the atomic medium. After a programmable delay 5t, a 
read pulse converts the state of atomic excitation into a 
field excitation, thereby generating one photon in a well- 
defined spatial and temporal mode 2. The quantum char- 




FIG. 1: Schematic of experiment for conditional generation 
of single photons. Write and read pulses sequentially prop- 
agate into a cloud of cold Cs atoms (MOT), generating the 
correlated output fields (1,2). A detection event for field 1 at 
D\A,iB leads to an approximate one-photon state for field 2, 
as confirmed with detectors D2A,2b- (P)BS - (polarization) 
beam splitter, SM - single-mode. The inset illustrates the 
relevant atomic level scheme. 

actcr of the (1, 2) fields is demonstrated by the observed 
violation of a Cauchy-Schwarz inequality for the ratio 
R of cross correlations to auto-correlations [111, narnely 
i? = 53±2 where i? < 1 for any classical field [20. 12lll2^. 

This greatly improved nonclassical correlation for pho- 
ton pairs for the (1, 2) fields enables the conditional gen- 
eration of single photons as in Refs. [I^, [H [3 El [H 
[T^ , but now with the photon 2 stored as an excitation in 
the atomic ensemble (15j . Given a first photon 1 from the 
write pulse, we trigger the emission of a second photon 
2 with the read pulse. To demonstrate the single-photon 
character of the field 2, we measure the three-fold cor- 
relation function 1^(12,12111) for detection of two pho- 
tons (I2, 12) from field 2 given a detection event li from 
field 1, where 1x1(12,12111) = 1 for coherent states and 
w(l2,l2|li) > 1 for any classical field. Experimentally, 
we find w(l2, 12II1) = 0.34 ± 0.06 for 5t = 150 ns, while 
w(l2, 12II1) = 0.24 ± 0.05 for 5t = 0, thereby taking an 
important step toward the creation of ideal single pho- 
tons for which w{\2, 12II1) 0. 

Figure ^ provides an overview of our experiment for 
producing correlated photons from an optically thick 
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sample of four-level atoms in a magneto-optical trap 
(MOT) m m. The ground states {|a); corre- 
spond to the 6Si/2,F — {4; 3} levels in atomic Cs, 
while the excited states {|e); |e')} denote the {6P3/2, -F" — 
4;6Fi/2,F = 4} levels of the D2,Di lines at {852; 894} 
nm, respectively. We start the protocol for single photon 
generation by shutting off all light responsible for trap- 
ping and cooling for 1/is, with the trapping light turned 
off approximately 300 ns before the re-pumping light in 
order to empty the F = 3 hyperfine level in the Cs 65i /2 
ground state, thus preparing the atoms in \a). During 

the "dark" period, the j'^ trial is initiated at time i^^-* 
when a rectangular pulse of laser light from the write 
beam, 150 ns in duration (FWHM) and tuned 10 MHz 
below the \a) — > |e) transition, induces spontaneous Ra- 
man scattering to level \b) via \a) — > |e) \b). The 
write pulse is sufficiently weak so that the probability to 
scatter one Raman photon into a forward propagating 

wavepacket ip^^\r,t'^-^'') is less than unity for each pulse. 
Detection of one photon from field 1 results in a "spin" 
excitation to level with this excitation distributed in 
a symmetrized, coherent manner throughout the sample 
of TV atoms illuminated by the write beam. 

Given this initial detection, the stored atomic excita- 
tion can be converted into one quantum of light at a user 
controlled time i^^'' = t'j^^ + St. To implement this con- 
version, a rectangular pulse from the read beam, 120 ns 
in duration (FWHM) and resonant with the \b) — > |e') 
transition, illuminates the atomic sample. This pulse af- 
fects the transfer \b) \e') \a) with the accompany- 
ing emission of a second Raman photon 2 on the |e') — > 

\a) transition described by the wavepacket "0*^^^ (r, tj-^"*), 

where the spatial and temporal structure of i/)^^'^) (r, t) 
are discussed in more detail in Ref. |2^. The trapping 
and re-pumping light for the MOT are then turned back 
on to prepare the atoms for the next trial j -I- 1 , with the 
whole process repeated at 250 kHz. 

The forward-scattered Raman light from the write, 
read pulses is directed to two sets of single-photon de- 
tectors {DiA,iB for field 1 and D2A,2B for field 2) [2^. 
Light from the {write, read) pulses is strongly attenu- 
ated (by ~ 10^) by the filters shown in Fig. ^ while 
the associated (1,2) photons from Raman scattering are 
transmitted with high efficiency (~ 80%) |20| . Detection 
events from Dia,ib within the intervals [t[^\t''.p + T] 
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and from D2A,2B within [t)j ,tj +T] are time stamped 
(with a resolution of 2 ns) and stored for later analysis. 
T — 200 ns for all of our measurements. 

For a particular set of operating conditions, we deter- 
mine the single pi and joint pi^rn 

event probabilities from 
the record of detection events at Dia,ib, D2A,2B, where 
{l,m) = 1 or 2. For example, the total singles probabil- 
ity pi for events at Dia , -Dis due to field 1 is found from 
the total number of detection events uia + riis recorded 
by DiA, DiB during the intervals [t^pit^p over Mtot 
repeated trials {j}, with then pi = {niA-^niB) / Mtot- To 
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FIG. 2: Left column (a) - (c) Normalized intensity correla- 
tion functions 31,1,^2,2,51,2 versus observed detection prob- 
abilities pi,p2, y'pip2, respectively. Right column (d) - (f) 
91,1, 92,2, gi,2 for joint detection versus gi, 92, ^91 92 for single 
detection, with qi,qi^rn referenced to the output of the MOT. 
Statistical uncertainties are indicated by the error bars. The 
full curves are from the model calculation described in the 
text with (ki,K2) = (0.17,0.90) and {\vit\^ , \v2b\^) = 0.006. 



determine for joint detections at Dia, Dib, we count 
the total number of coincidences Nia^ib recorded by 
DiA, DiB, with thenpi_i = NiA,iB/Mtot- P2,2 is found in 
an analogous fashion using events from D2A, D2b- Joint 
detections between the (1, 2) fields are described by pi,2, 
which is determined by summing coincidence events be- 
tween the four pairs of detectors for the (1, 2) fields (e.g., 
between pairs Dia, D2a)- 

From (j>i,pi^m) we derive estimates of the normalized 
intensity correlation functions gi^rm where gi^m = 1 for 
coherent states. For example, the auto-correlation func- 
tion gi i = pi.i/ (piaPib) for field 1, and similarly for the 
functions 52,2,51,2 for the auto-correlation of field 2 and 
the cross correlation between fields (1,2). The first col- 
umn in Figure 121 displays 51,1,52,2, and 51,2 as functions 
oipi,P2, and ^JpiP2 A virtue of 5;,™ is its indepen- 
dence from the propagation and detection efficiencies. In 
the ideal case, the state for the fields (1.2) is |fiL l2l[2i | 

1*12) - |00) + Vxlll) + x|22) + 0(x^/') , (1) 

where ,Jx is the excitation amplitude for field 1 in each 
trial of the experiment. For x ^ 1, 5i,i = 52,2 — 2 
and 51,2 = 1-1- 1/x- By contrast, for reasons that we 
will shortly address, our measurements in Fig. [3 give 
51,1 ~ 1.7 and 52,2 — 1-3, with 51,2 exhibiting a sharp rise 
with decreasing y^PiP2, but with considerable scatter. 

To provide a characterization of the field generation 
that is independent of the efficiency of our particular 
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detection setup, we convert the photodetection proba- 
bilities ipi,pi^m) to the quantities {qi,qi^m) for the field 
mode collected by our imaging system at the output of 
the MOT. Explicitly, for single events for fields (1,2), 
we define qi = pi/ai, q2 = V2I d^-, while for joint events, 
<7i,i = V\.\la\,q2.2 = 2/027 91,2 = Pi,2/ckia2, where ai 
give s the collection, propagation, and detection efficiency 
[2g. The second column in Fig. |21 displays the measured 
dependence of qi^m for joint events versus q\,qi, ^q\qi 
for single events over a range of operating conditions. As 
expected from Eq. q\;\_,qi^i exhibit an approximately 
quadratic dependence on qi,q2, while qi^2 would be linear 
for y/qiq2 ^ 1 in the ideal case. 

In our experiment there are a number of imperfections 
that lead to deviations from the ideal case expressed by 
1^12) [EHEI^' To capture the essential aspects, we have 
developed a simple model that assumes the total fields 
(1,2) at the output of the MOT consist of contributions 
from I $12), together with background fields in coherent 
states |wi^2)- Operationally, increases inpi,p2 sue accom- 
plished by way of increases in the intensity of the write 
beam, with only minor adjustments to the read beam. 
Hence, we parameterize our model by taking x — ktuPi 
with the (scaled) amplitude of the write beam. 

Since important sources of noise are light scattering from 
the write and read beams and background fiuorescence 
from uncorrelated atoms in the sample 12J|, we assume 
that vi^2 = y/Ki^2Vw We further allow for fixed incoher- 
ent backgrounds vu, V2b to account for processes that do 
not depend upon increases in the write intensity. 

With this model, it is straightforward to compute the 
quantities that appear in Figs. EH The parameters 
{ki,K2) = (0.17, 0.90) and {\vit\^, \v2b\^) = 0.006 are ob- 
tained directly by optimizing the comparison between the 
model results and our measurements of normalized cor- 
relation functions (e.g., gi^i vs. 51,2) without requiring 
absolute efficiencies, ki = 0.17 implies that the photon 
number for "good" events associated with |$i2) exceeds 
that for "bad" (background) events from \vi) by roughly 
6— fold for detection at Dia, Dib- For the curves in Fig. 
121 we must also obtain the efficiencies Pi, rji that convert 
expectation values for normally ordered photon number 
operators hi for fields / = (1,2) in the model into the 
various {puPi^m) and {qi,qLrn) (e.g., pi = Pi{hi),qi = 
Vii^i), 11.2 ^mmi- nih2- ))■ Ideally /3i = ai and 77; = 1; 
we find instead {(3i,r]i) = (0.013,0.15), where we take 
/3i = f32 and rji — r]2 for simplicity. Among various can- 
didates under investigation, values Pi < ai,r]i < 1 can 
arise from inherent mode mismatching for capturing col- 
lective emission from the atomic ensemble |2J|. 

Independent of the absolute efficiencies, we can uti- 
lize the results from Fig. |21to address directly the ques- 
tion of the nonclassical character of the (1, 2) fields by 
following the pioneering work of Clauser lllj. The cor- 
relation functions gi^m for fields for which the Glauber- 
Sudarshan phase-space function ip is well-behaved (i.e., 
classical fields) are constrained by the inequality R = 
[31,2] Vffi, 152,2 < 1 [Hill- In Fig. El we plot the ex- 




FIG. 3: Ratio R= [§1,2] /gi,ig22 versus the normalized cross 
correlation §1,2, where R > 1 for manifestly quantum (non- 
classical) fields. The points are from our experiment with 
statistical uncertainties indicated by the error bars. The 
full curve is from the model calculation with (ki,K2) and 
i\vitf,\v2b\^) as in Fig.H 



perimentally derived values for i? as a function of the 
degree of cross-correlation 31^2 jlsj. As compared to pre- 
vious measurements for which R = 1.84 ± 0.06 and 
R = 1.34 ± 0.05 fU, we have now achieved i? > 1, with 
i? = 53 ± 2 for the largest value of 51.2. In Figs. [21and|3| 
as well as 01 to follow, all points are taken with St — 150 
ns, except the points at gi,2 — 10, which have 6t — 0. 

This large degree of quantum correlation between the 
(1,2) fields suggests the possibility of producing a sin- 
gle photon for field 2 by conditional detection of field 1. 
To investigate this possibility, we consider the three-fold 
correlation function i(;(l2, 12II1) for detection with the 
setup shown in Fig. ^ namely 



W(l2,l2|ll) 



^ P<^Hl2,l2|ll) 



(2) 



where p^'^^ (I2, 12II1) is the conditional probability for de- 
tection of two photons (12,12) from field 2 conditioned 
upon the detection of an initial photon li for field 1, and 
p'-'^Hlajli) is the probability for detection of one photon 
I2 given a detection event li for field 1. Bayes' theorem 
allows the conditional probabilities in Eq. |21to be written 
in terms of single and joint probabilities p*-*^^ for /c— fold 
detection, so that 

pW(li)p(3)(i,, 

"^'^''^''^^= b(^)(li,l2)P • 

Fields with a positive-definite ip must satisfy the Cauchy- 
Schwarz inequality w(l2, 12 1 li) > 1- Indeed, for indepen- 
dent coherent states, w ~ 1, while for thermal beams, 
w ^ 2. By contrast, for the state |<I>i2) of Eq. ^ 
w = 4% <C 1 for small x, approaching the ideal case 
w — *■ for a "twin" Fock state |lil2)- 

From the record of photo-detection events at 
DiA.iB, D2A.2B^ we calculate estimates of the various 
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for field 2). In this case, the (1,2) fields should be sta- 



FIG. 4: Three-fold correlation function w(l2,l2|li) for de- 
tection event li for field 1 followed by two events (I2, 12) 
for field 2 versus the normalized cross correlation gi,2. (a) 
Wi,j{l2, l2|f 1) for events (li)i and (I2, f2)j from different tri- 
als i ^ j together with points for §2.2- = 52,2 for statis- 
tically independent trials, (b) 'u;i_i(f2, f2|f 1) for events from 
the same trial i. Wi^i < 1 for sub-Poissonian fields in support 
of the single-photon character of field 2. Statistical uncertain- 
ties are indicated by the error bars. The full curves are from 
the model calculation with (ki,K2) and {\vit\^ ,\v2b\^) as in 
Figs. HEl 



probabilities appearing in Eq. 13 with the results of this 
analysis shown in Fig. 01 Part (a) examines the quantity 
U'ij (l2, 12II1) obtained from events taken from different 
trials i j for the (1,2) fields (i.e., detection li in trial 
i for field 1 followed by two detections (I2, 12) in trial j 



tistically independent [23, so that ^1,^(12,12111) = 92,2- 
Hence, we also superimpose 32,2 from Fig. |21 and find 
reasonable correspondence within the statistical uncer- 
tainties (in particular, (I2, Igjli ) ^ 1), thereby val- 
idating our analysis techniques |25j| . No corrections for 
dark counts or other backgrounds have been applied to 
the data in Fig. El(nor indeed to Figs. [21I3I)- 

Fig. E (b) displays Wi,i(l2, 12II1) for events from the 
same experimental trial i for the (1,2) fields. Signifi- 
cantly, as the degree of cross-correlation expressed by 
.gi,2 increases (i.e., decreasing x), Wi,i(l2, 12II1) drops 
below the classical level of unity, indicative of the sub- 
Poissonian character of the conditional state of field 
2. For St = 150 ns, w^,»(l2, 12II1) = 0.34 ± 0.06 for 
gi,2 = 7.3, while with 5t = 0, w^,,(l2, 12II1) = 0.24±0.05 
for gi.2 = 10.2. Beyond the comparison to our model 
shown the figure, empirically we find that Wi,i(l2, 12II1) 
is well approximated by 51.152,2/51,2, as in the ideal case 
of Eq. n However, independent of such comparisons, we 
stress that the observations reported in Fig. ^represent 
a sizable nonclassical efi'ect in support of the conditional 
generation of single photons for field 2. 

In conclusion, our experiment represents an important 
step in the creation of an efficient source of single photons 
stored within an atomic ensemble, and thereby towards 
enabling diverse protocols in quantum information sci- 
ence 0,0, IE Q' Our model supports the hypothesis 
that the inherent limiting behavior of ^^,^(12, 12II1) be- 
low unity is set by the efficiency 77;, which leads to pro- 
hibitively long times for data acquisition for x ^ 0.04, 
corresponding to the smallest value of Wi^i in Fig. 0] We 
are pursuing improvements to push rji ~ 0.15 — > 1, in- 
cluding in the intrinsic collection efficiency following the 
analysis of Ref. '24'|. Dephasing due to Larmor preces- 
sion in the quadrupole field of the MOT limits St < 300 
ns, which can be extended to several seconds in optical 
dipole or magnetic traps [23|. 
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